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FOREWORD 

This  document  presents  the  results  of  an  experimental  and  analytical 
investigation  into  the  characteristics  of  the  three  dimensional  shock 
wave  turbulent  boundary  layer  interaction.  The  study  was  conducted  by 
the  High  Speed  Aero  Performance  Branch  (FXG),  Aeromechanics  Division, 

Air  Force  Flight  Dynamics  Laboratory,  Wright  Patterson  Air  Force  Base, 

Ohio.  The  work  concludes  an  in-house  and  contracted  research  program 
and  was  performed  under  Project  13G6  "Aeroperformance  and  Aeroheating 
Technology",  Task  136603  "Aerodynamic  Heating  to  Military  Vehicles". 

This  report  covers  experimental  work  conducted  from  January  1974  to 
September  1976  and  concludes  work  unit  13660329  "Design  Techniques  to 
Predict  and  Minimize  Interference  Heating’. 

The  contracted  research  effort  was  as  follows.  Contract  F33615-73- 
C-3133  was  with  Princeton  University,  Department  of  Aerospace  and  Mechanical 
Sciences,  Princeton,  New  Jersey  and  covered  the  period  of  April  1973 
through  November  1974.  This  contractor  investigated  and  reported  on  heat 
transfer,  surface  pressure,  and  flow  field  data  for  Mach  2.95.  Contract 
F33615-73-C-3046  with  McDonnell  Douglas  Corporation,  St.  Louis,  Missouri 
covered  the  period  of  January  1973  through  April  1974.  This  contractor 
investigated  and  reported  on  pressure  data  and  Heat  Transfer  at  Mach  3.71. 
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SUMM.ARY 

The  characteristics  of  the  three  dimensional  shock  wave  turbulent 
boundary  layer  interaction  were  investigated  through  an  extensive 
experimental  test  program.  This  program  encompassed  Mach  numbers  fran 
2.95  to  5.85,  Reynolds  numbers  of  1.5  x 10^  to  28  x 10^  per  foot,  and 
boundary  layer  thicknesses  a^'  the  fin  leading  edge  from  0.13  to  6 inches. 
Fach  test  collected  heat  transfer,  pressure,  and  oil  flow  data  at  selected 
fin  deflection  angles  ranging  from  0°  to  20°.  The  objective  was  to 
generate  a data  base  covering  a wide  enough  range  of  parameters  to  allow 
the  formulation  of  a set  of  inclusive  empirical  prediction  techniques  for 
the  major  characteristics  of  the  interaction. 

This  report  presents  the  results  of  the  experimental  program  and  the 
correlations  developed  from  the  data.  It  was  found  that  most  of  the 
interaction  characteristics  were  dependent  on  the  value  of  the  shock 
strength  (M^  Sin  0)  and  on  the  distance  aft  of  the  fin  leading  edge 
measured  in  boundary  la^er  thicknesses  (X/6).  The  correlations  developed 
for  the  peak  pressure  and  peak  heating  were 

= (M^  Sin  e)  P and  ~ = n..  (M  Sin  6 - 1)  + 0.75 

Tpp 

where  Op  and  are  functions  of  X/6.  At  large  X/6  it  was  observed  that 
Op  = 2.40  and  n^^  = 4.75.  Both  represent  values  substantially  higher  than 
2-D  peak  values.  By  correlating  2-0  and  3-D  plateau  pressures  in  the 
separated  region  Spaid's  conclusion  that  the  separation  phenomena  is 
independent  of  the  method  used  to  generate  it  was  substantiated  and 
extended  to  3-D  interactions.  The  extent  of  the  interaction  region  was 
determined  by  correlating  the  location  of  the  separation  line.  It  was 
observed  that  the  separation  line  could  be  described  by  a set  of  hyperbolic 
curves  when  the  correct  coordinate  system  was  established.  These  hyperbolic 
curves  were  defined  by  correlating  the  slope  and  interceots  of  their 
asymptotes  as  a function  of  Sin  0. 


X 


AFFDL-TR-77-10 


Complementing  the  above  results  are  correlations  on  the  shock  and 
separation  pressure  ratios,  the  location  of  the  peak  pressure  and  heat 
transfer  rate,  and  the  pressure  distribution  through  the  interaction 
region. 
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SECTION  I 
INTRODUCTION 


Sustained  flight  at  high  Mach  numbers  presents  a variety  of  problems 
to  the  engineer.  A fundamental  and  highly  influencial  phenomenon  is  that 
of  the  shock  wave  boundary  layer  interaction.  This  flow  field  interaction 
is  characterized  by  boundary  layer  separation  and  by  localized  areas  of 
high  pressure  and  high  aerodynamic  heating  which  impact  on  the  design  of 
empennages  and  control  surfaces. 

The  two  dimensional  interaction  has  received  much  experimental 
attention  and  commonly  occurs  on  flight  vehicles  at  control  surface 
hinge  lines.  The  three  dimensional  interaction,  characteristically 
generated  by  deflected  fins  and  boundary  layer  diverters,  has  received 
less  attention  although  it  produces  more  severe  pressure  and  heating 
problems.  Moreover,  the  published  data  on  3-D  interactions  generally 
presents  conflicting  trends  in  peak  pressure  and  peak  heating  correlations. 
This  undoubtedly  is  a result  of  the  complexity  of  the  interaction  and  the 
fact  that  no  single  investigator  has  obtained  data  over  a sufficient  range 
of  Mach  numbers,  Reynolds  numbers, and  fin  deflection  angles. 

The  Air  Force  Flight  Dynamics  Laboratory  (AFFDL)  has  an  ongoing 
program  involved  in  the  collection  and  generation  of  3-D  interaction 
data.  The  objective  is  to  obtain  data  over  a wide  enough  range  of  the 
above  parameters  to  allow  the  formulation  of  an  inclusive  prediction 
technique  for  the  major  interaction  characteristics.  This  report  presents 
the  current  results  of  the  program. 

Table  1 presents  the  data  base  on  which  the  correlations  in  this 
report  were  made.  Oskam,  Bogdonoff,  and  Vas  (Reference  1)  carried  out 
a study  at  Mach  2.95  in  the  Princeton  University  high  Reynolds  number 
tunnel.  In  this  investigation  two  model  configurations  were  employed 
to  study  the  effect  of  boundary  layer  thickness  on  the  interaction.  A 
boundary  layer  thickness  of  0.13  inch  at  the  fin  leading  edge  was  obtained 
by  mounting  a fin  on  a sharp  flat  plate  suspended  between  the  tunnel  side 


1 


TABLE  1 
DATA  BASE 


AFFDL-TR-77-10 


walls.  The  second  configuration  mounted  the  fin  directly  on  the  floor 
of  the  tunnel  to  obtain  a boundary  layer  thickness  of  0.55  inch  at  the 
fin  leading  edge.  A similar  study  was  made  in  the  NASA  Langley  Unitary 
Plan  Tunnel.  Token  (Reference  2)  carried  out  a study  at  Mach  3.71  in 
which  a fin  15  inches  high  and  32  inches  long  was  mounted  on  the  tunnel 
side  wall  in  a 6 inch  boundary  layer.  In  a follow  up  study  (Reference  3) 
AFFDL  mounted  the  same  fin  on  a sharp  flat  plate  6 feet  long,  and 
suspended  the  plate  from  the  tunnel  side  wall.  Data  were  obtained  at 
Mach  numbers  of  3.0  and  3.7  with  boundary  layer  thicknesses  ranging  from 

0. 37  to  0.43  inches. 

A large  amount  of  data  have  been  obtained  in  the  AEDC-VKF  Tunnels  A 
and  B.  A small  flat  plate  model  18  inches  long  and  10  inches  wide  was 
employed  in  Tunnel  A.  A fin  7.5  inches  long  and  3 inches  high  was  mounted 
on  the  plate.  AFFDL  obtained  data  (Reference  4)  at  Mach  numbers  of  3.0, 
3.75,  and  4.5  with  this  configuration.  The  same  fin  was  then  mounted  on  a 
plate  24  inches  square  and  installed  in  Tunnel  B.  The  author  obtained 
data  on  this  model  at  Mach  4.75  and  5.04.  This  previously  unpublished 
data  was  presented  at  the  American  Institute  of  Aeronautics  and  Astro- 
nautics (AIAA)  15th  Aerospace  Sciences  Meeting  (Reference  5). 

The  model  employed  in  the  Tunnel  A test  was  also  installed  in  the 
AFFDL  Mach  6 High  Reynolds  Number  Facility.  Data  (Reference  6)  was 
obtained  in  this  tunnel  at  a nominal  Mach  number  of  5.85. 

The  range  in  Mach  number  and  Reynolds  number  which  the  data  base 
covers  has  made  it  possible  to  correlate  sets  of  data  which  previously 
were  considered  to  present  conflicting  data  trends. 

1.  BASIC  3-D  INTERACTION  CHARACTERISTICS 

The  typical  3-0  interaction  generator  is  shown  in  Figure  1.  It 
consists  of  a sharp  flat  plate  instrumented  with  thermocouples  or  pressure 
taps,  and  a fin  mounted  normal  to  the  plate  which  may  be  set  to  various 
deflection  angles.  Also  shown  in  Figure  1 are  typical  pressure  and  heat 
transfer  profiles  taken  in  a plane  normal  to  the  undisturbed  flow  direction. 
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The  pressure  distribution  may  be  divided  into  two  distinct  regions 
by  the  fin  oblique  shock  wave.  The  outer  region  extends  from  onset  of 
the  interaction  region  to  the  shock  wave  and  is,  at  sufficient  shock 
strength,  a region  of  separated  flow.  The  flow  in  this  region  has  distinct 
two  dimensional  characteristics.  As  will  be  shown  later  the  separation 
and  plateau  pressure  ratios  may  be  predicted  with  2-D  data  correlations. 

The  inner  region,  extending  from  the  shock  wave  to  the  fin  surface,  is 
characterized  by  the  sharp  peak  in  pressure  which  lies  close  to  the  fin. 

At  sufficient  shock  strength  the  flow  in  this  region  is  dominated  by  an 
imbedded  vortex  which  produces  a trough  or  low  pressure  region  under  the 
shock  wave  and  peak  pressures  in  excess  of  oblique  shock  predictions. 

The  heat  transfer  distribution  is  characterized  by  a small  peak  at 
the  line  of  boundary  layer  separation,  followed  by  a smooth  rise  to  the 
shock  wave  location.  At  the  shock  location  the  slope  steepens  abruptly 
as  the  fin  is  approached.  The  large  inner  peak  occurs  at  approximately 
the  same  location  as  the  pressure  peak. 

When  sufficient  shock  strength  is  attained  to  produce  separation 
the  oil  flow  data  have  characteristics  as  shown  in  Figure  1.  The  peak 
pressure  and  heating  occur  along  the  line  of  streamline  divergence  near 
the  fin.  Maximum  streamline  deflection  occurs  at  thp  shock  location  as 
determined  by  oblique  shock  relations.  The  line  of  boundary  layer 
separation  is  observed  to  be  the  inner  edge  of  the  oil  accumulation  line 
and  onset  is  the  point  at  which  the  undisturbed  streamlines  first  begin 
to  curve. 
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SECTION  II 
PRESSURE  PROFILES 

The  features  of  the  pressure  profile  introduced  in  Section  I-l  which 
will  be  discussed  are  the  peak,  shock,  plateau,  and  separation  pressure 
magnitudes  and  their  locations  in  the  interaction  region.  The  correlation 
parameters  found  to  be  most  successful  are  the  shock  strength,  Sin  6, 
and  the  nondimensional  distance,  "X/6".  The  parameter  Sin  6 is  the 
component  of  the  free  stream  Mach  number  normal  to  the  fin  oblique  shock 
wave,  and  is  calculated  from  oblique  shock  relations.  This  parameter  is 
useful  as  a measure  of  the  strength  of  the  interaction  and  as  a trans- 
formation variable  in  applying  2-D  correlations  to  the  3-D  data.  The 
parameter  "X/6"  is  the  ratio  of  the  distance  downstream  of  the  fin  leading 
edge  to  the  boundary  layer  thickness  at  the  fin  leading  edge.  The  inter- 
action region  grows  in  strength  and  extent  with  distance  aft  of  the  fin 
leading  edge.  The  pressure  in  the  plateau  region  remains  nearly  constant 
with  X/6,  however  the  peak  pressure  increases  until  it  is  greater  than 
oblique  shock  predictions. 

1 . PEAK  PRESSURE  RATIO 

The  peak  pressure  ratios  were  investigated  by  plotting  on  log-log 
paper  the  peak  pressures  as  a function  of  Sin  6.  At  each  Mach  number 
and  for  each  X position  the  result  is  a straight  line  passing  through 
PpK/Ppp  = 1 and  M^  Sin  6 = 1.  The  slope  of  the  line  seemed  to  vary  with 
Mach  number  and  X,  however  it  wa<;  found  that  the  apparent  Mach  number 
dependence  could  be  removed  by  considering  the  dependence  on  X/6  rather 
than  X.  Figure  2 shows  several  representative  plots  of  the  peak  pressure 
data  as  a function  of  M^  Sin  6 and  X/6.  They  indicate  a relationship  of 
the  form 


{M  Sin  e) 

' oo  ' 


(1) 


where  Op  is  the  slope  of  the  curves  in  Figure  2 and  is  a function  of 

X/6., 
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All  of  the  available  peak  pressure  data  were  plotted  in  the  format 

of  Figure  2 and  the  value  of  n„  calculated.  The  function  n„  derived  is 

P P 

shown  in  Figure  3.  Figure  3 shows  that  at  small  X/6  the  peak  pressure 
is  lower  than  oblique  shock  predictions  regardless  of  shock  strength. 
However,  the  oblique  shock  values  are  attained  within  ten  boundary  layer 
thicknesses  of  the  fin  apex.  At  X/6  greater  than  ten  the  oblique  shock 
values  are  exceeded.  It  should  also  be  noted  that  the  value  of  Op 
approaches  a constant  value  of  about  2.4  at  large  X/6.  This  shows  that 
although  the  pressure  continues  to  increase  with  shock  strength  at  each 
X/6,  the  rate  at  which  this  rise  occurs  approaches  a fixed  value.  At 
X/6  > 30  the  peak  pressure  function  then  is 

(M„  Sin  6)^-'^  (2) 

^FP 

for  any  X/6. 

The  peak  pressure  occurs  along  a line  which  mai  be  closely  approximated 
by  a ray  from  the  fin  leading  edge.  The  angle  between  this  ray  and  the 
free  stream  direction  has  been  correlated  in  Figure  4 with  the  value 
predicted  by  Token's  expression  (Reference  2). 

= 0.24  (6-a)  + a (3) 

2.  PRESSURE  RATIO  UNDER  THE  SHOCK 

Figure  5 presents  the  development  of  the  pressure  distribution  with 
shock  strength.  At  low  shock  strength,  below  that  required  for  incipient 
separation,  the  pressure  rises  smoothly  from  onset  to  the  peak  value. 

When  incipient  separation  is  reached  a plateau  region  begins  to  form 
outside  the  shock  wave  location.  Increased  shock  strength  causes  the 
formation  of  a trough  under  the  shock  wave  location.  The  pressure  ratio 
in  this  trough  varies  with  shock  strength  and  distance  aft  of  the  fin 
leading  edge. 
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Figure  3.  Exponent  in  Peak  Pressure  Correlation 
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The  first  attempt  to  correlate  the  shock  pressure  in  the  trough  was 
based  on  the  oil  flow  data  which  shows  that  the  boundary  layer  flow  inside 
the  shock  wave  moves  outward  under  the  shock  and  into  the  separated 
region.  It  was  thought  that  the  trough  pressure  level  might  be  determined 
by  the  pressure  rise  across  the  shock  external  to  the  boundary  layer.  It 
was  also  observed  that  at  large  X/5  the  pressure  level  in  the  trough  varied 
with  shock  strength  in  a manner  similar  to  that  of  the  unit  Reynolds 
number  ratio  across  the  shock.  That  is,  the  pressure  rises  with  Sin  6 
to  some  maximum  value  and  then  drops  off  again  as  Sin  0 increases. 

The  maximum  pressure  occurs  at  the  same  shock  strength  as  the  maximum 
Reynolds  number  ratio.  However  no  conclusive  correlation  could  be  found 
between  the  shock  pressure  level  and  either  the  Reynolds  number  or  the 
pressure  differential  across  the  shock. 

A second  and  more  successful  attempt  to  correlate  the  trough  pressure 
was  based  on  Token's  (Reference  2)  work.  Token  suggested  that  the  peak 
heating  rate  could  be  correlated  with  stagnation  point  theory  and  developed 
equations  for  the  heating  distribution  between  the  peak  and  shock  wave 
locations.  The  equations  Included  a pressure  gradient  term  involving  the 
peak  and  shock  pressure  ratios.  The  peak  Stanton  number  was  therefore 
plotted  as  a function  of  (Ppj^  - P5^)/Ppp  shown  in  Figure  6.  The  function 
shows  a dependence  on  X/6,  however  the  limited  amount  of  data  prevents 
firm  conclusions  on  it's  effects  from  being  drawn.  Each  point  on  this 
figure  requires  accurate  heat  transfer  and  pressure  data  at  the  peak  and 
shock  locations  and  that  the  shock  strength  be  large  enough  to  produce 
the  trough.  A small  percentage  of  the  data  base  satisfied  all  these 
requirements.  The  data  in  Figure  6 indicate  that  this  correlation 
together  with  the  peak  pressure  and  peak  heating  correlations  can  be  used 
to  determine  the  shock  pressure,  however  more  data  is  required  to  increase 
its  utility, 

3.  PLATEAU  PRESSURE  RATIO 

The  flow  properties  between  onset  and  the  shock  wave  in  the  3-D 
interaction  are  similar  to  those  in  the  separated  region  of  a 2-D  interactions 
If  Sin  0 is  used  as  the  Mach  number  for  the  interaction  then  the  plateau 
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and  separation  pressures  can  be  predicted  by  2-D  data  correlations. 

Figure  7 shows  the  2-D  plateau  pressure  correlation  curves  and  2-0  and 
3-D  data  superimposed.  The  plateau  pressure  is  independent  of  X,  6,  and 
Reynolds  number  and  depends  only  on  the  shock  strength. 

The  plateau  pressure  is  a result  of  the  boundary  layer  separation 
process  and  its  magnitude  is  independent  of  the  method  used  to  generate 
the  separation.  Data  from  2-D  shock  impingements,  ramps,  3-D  fin  shock 
impingement,  steps,  and  jet  injection  all  correlate  well  in  the  plateau 
region. 

4.  SEPARATION  PRESSURE  RATIO 

The  boundary  layer  separation  line  in  the  3-D  interaction  can  be 
identified  in  oil  flow  photographs  as  the  inner  edge  of  the  oil  accumulation 
line  produced  by  the  convergence  of  the  freestream  surface  streamlines  as 
shown  in  Figure  1.  The  surface  pressure  at  this  location  agrees  well  with 
2-0  separation  theory  which  predicts  separation  at  Tit  of  the  plateau 
pressure  rise.  Figure  8 presents  the  separation  pressure  data  superimposed 
on  a curve  representing  lit  of  the  plateau  correlation  curve. 

Accurately  locating  the  separation  line  without  oil  flow  dat.i  is 
difficult.  It  was  found  that  the  separation  pattern  around  half  cones  on 
a flat  plate  formed  hyperbolic  curves  (Reference  12).  This  approach  was 
applied  to  the  fin  data.  A coordinate  system  was  set  up  with  the  fin 
shock  wave  on  the  X axis  and  the  fin  leading  edge  at  a distance  "a"  from 
the  Y axis.  In  this  coordinate  system  the  separation  lines  form  hyperbolic 
curves.  As  shown  in  Figure  9,  a hyperbolic  curve  is  defined  by  its  asymptotes. 
The  asymptotes  may  be  described  by  their  slope  "b/a"  and  either  "a"  or  "b"; 
therefore,  the  separation  line  may  be  defined  by  two  numbers  which  must  be 
correlated  in  order  to  predict  separation  location.  The  slope  of  the 
asymptotes  correlates  well  with  Sin  9 as  shown  in  Figure  10.  This 
figure  shows  that  the  slope  ^ approaches  zero  at  Sin  0 1.3.  At  this 

shock  strength  the  shock  wave  and  the  separation  line  are  coincident.  In 
reality  separation  no  longer  exists.  The  pressure  plateau  also  ceases  to 
exist  and  the  peak  pressure  approaches  the  oblique  shock  value.  For  shock 
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Figure  7.  Plateau  Pressure  Correlation 
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strengths  less  than  1.3  the  pressure  distribution  is  characterized  by  a 
smooth  rise  from  the  freestream  value  to  the  peak  or  oblique  shock  value. 
The  value  of  either  "a"  or  "b"  are  more  difficult  to  correlate.  Because 
the  slope  of  the  asymptotes  is  so  small  the  variation  of  "a"  is  large  with 
small  changes  in  b/a.  Figure  11  shows  the  correlation  obtained  for  "a" 
as  a function  of  Sin  9.  As  Sin  0 approaches  1.3  "b/a"  approaches 
zero  and  "a"  approaches  infinity.  This  infers  that  the  condition  for 
incipient  separation  is  that  Sin  0 equal  1.3.  Korkegi's  (Reference  13) 
incipient  separation  condition  was 


M = 0.3 

» inc. 


(4) 


If  the  incipient  separation  condition  Sin  0 = 1.3  is  transformed  to 
Korkegi's  form  the  results  are 


for  3 < M <6 
— 00  — 


0*423  < <0.438 


(5) 


Equation  5 then  predicts  separation  at  a slightly  higher  fin  deflection 
then  Equation  4. 


Figure  12  presents  the  separation  coordinates  as  a function  of  Sin  9. 
These  curves  were  generated  using  curve  fits  of  Figures  10  and  11,  and  the 
equation 


(6) 


5.  ONSET  OF  PRESSURE  RISE 

The  point  of  onset  was  defined  in  Section  I-l  by  the  point  in  the 
oil  flow  data  at  which  the  undisturbed  surface  streamlines  first  begin 
to  curve.  This  point  may  be  located  on  the  pressure  profiles  by  fairing 
a line  through  the  pressure  rise  to  the  point  of  separation.  As  shown  in 
Figure  1 the  point  at  which  this  line  intersects  P/Ppp  = 1 is  taken  as 
the  point  of  onset.  Having  calculated  the  location  of  separation  as 
shown  in  Section  II-4  the  location  of  onset  may  be  calculated  by  correlating 
the  dP/dY  from  onset  to  separation.  The  magnitude  of  dP/dY  was  observed 
to  be  a function  of  shock  strength  and  distance  aft  of  the  fin  leading 
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edge.  As  shown  in  Figures  13  to  17  AP/AY  was  plotted  as  a function  of 
Sin  0 for  each  Mach  number  where  sufficient  data  were  available.  AP/AY 
is  found  to  be  a linear  function  of  shock  strength  at  each  station  within 
the  interaction  region.  The  slope  of  this  function  is  dependent  on  distance 
aft  of  the  fin  leading  edge  and  on  boundary  layer  thickness  as  shown  in 
Figure  18.  When  constructing  a pressure  profile  like  the  one  shown  in 
Figure  1 the  location  of  onset  is  determined  by  first  locatimi  separation 
and  then  extending  the  profile  to  onset  using  an  estiniation  of  AP/AY  from 
Figure  18. 

6.  PRESSURE  DISTRIBUTIONS 

Discussion  of  the  pressure  distributions  will  be  confined  to  data 
taken  at  shock  strengths  greater  than  Sin  6 = 1.3.  These  are  characterized 
by  the  vortex  dominated  flow  where  the  peak  pressure  is  greater  than  the 
oblique  shock  values,  a trough  is  observed  at  the  shock  location,  and  a 
distinct  plateau  region  exists.  Because  of  the  complexity  of  the  3-D 
interaction  the  pressure  distributions  must  be  broken  up  into  several 
sections  for  any  correlation  to  be  obtained. 


The  pressure  distribution  is  first  separated  Into  two  main  parts  by 
the  shock  wave.  The  inner  region  consists  of  a smooth  and  sharp  rise  from 
the  shock  pressure  to  the  peak  value  near  the  fin.  No  correlation  with 
2-D  data  has  been  found  here.  The  outer  region  however  has  distinct  2-D 
characteristics  as  described  in  Sections  II-3  and  II-4.  This  outer  region, 
extending  from  onset  to  the  shock  location,  is  divided  into  two  subregions 
by  the  point  of  boundary  layer  separation.  The  pressure  distribution 
corresponds  to  that  of  a 2-0  interaction  of  equal  shock  strength. 


a.  Inner  Region 


The 

inner  region 


normalized  coordinates  for  plotting  the  distribution  in  the 
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Figure  13.  Onset  Pressure  Gradient  at  Mach  2.95 


Figure  15.  Onset  Pressure  Gradient  at  Mach  3.75 


idient  at  Mach  5.04 


Figure  18.  Onset  Pressure  Gradient  Correlati 
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Figure  19  presents  the  correlation  for  the  inner  region  pressure  distribution. 
Data  scatter  is  obvious  and  to  be  expected  from  such  an  elementary  correlation 
method.  However  the  peak  pressure  and  its  location  can  be  accurately 
determined  from  Figures  3 and  4 and  Figure  19  is  then  adequate  for  producing 
a fast  and  simple  engineering  approximation  to  the  distribution  between 
the  peak  and  the  shock. 


b.  Outer  Region 

The  outer  region  which  extends  from  the  shock  wave  to  onset  is 
divided  into  two  subregions.  Subregion  1 extends  from  the  shock  to  the 
separation  line.  The  normalized  coordinates  employed  in  correlating  the 
data  in  this  subregion  are 


P 

P 


" ^FP 
PL  " *^FP 


and 


Y 


Figure  20  presents  the  pressure  distribution  in  this  subregion.  The  data 
correlation  is  good  only  to  (Y^,p  - Y)/Yj^p  = 0.32  which  may  be  used  as 
the  location  of  the  plateau  for  engineering  purposes.  Between  the  plateau 
and  the  shock,  the  pressure  distribution  is  a strong  function  of  shock 
strength.  At  Sin  9 = 1.3  the  pressure  remains  fairly  constant  between 
the  plateau  and  shock  locations.  As  Sin  6 increases  the  pressure 
begins  to  fall  off  into  a trough  at  the  shock  which  deepens  with  Sin  6 
and  X. 

Subregion  2 extends  from  the  separation  line  to  the  onset  of  the 
interaction  region.  The  normalized  coordinates  employed  in  the  data 
correlation  were 


and 


The  correlation  is  shown  in  Figure  21.  The  point  of  onset  may  be  measured 
from  oil  flow  data  at  the  point  where  the  freestream  surface  streamlines 
first  begins  to  curve  or  it  can  be  calculated  from  Figure  18. 
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Figure  19,  Pressure  Distribution,  Peak-to-Shock 
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The  pressure  distribution  has  been  correlated  from  onset  to  the  plateau 
and  from  the  shock  to  the  peak  location.  Work  is  still  required  to 
increase  the  utility  of  Figures  6 and  18  by  additional  data,  however  the 
characteristics  important  from  an  engineering  standpoint  have  been  defined. 
These  are  the  location  and  magnitude  of  the  pressure  at  separation,  at  the 
plateau  and  at  the  peak. 
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SECTION  III 

HEAT  TRANSFER  PROFILES 


A typical  heat  transfer  profile  was  introduced  in  Section  I-l.  The 
dominant  characteristics  are  a small  peak  at  the  line  of  separation 
followed  by  a nearly  linear  rise  to  the  shock  location.  At  the  shock  a 
sharp  change  in  slope  occurs  and  the  heating  rates  continue  to  rise  to 
the  large  inner  peak. 


1 . PEAK  HEAT  TRANSFER 

It  was  shown  in  Section  II-l  that  the  peak  pressure  data  taken  in 
flows  having  large  differences  in  boundary  layer  thickness  could  be 
correlated  through  the  term  "X/6".  This  type  of  correlation  has  been 
observed  to  hold  for  the  peak  heating  data  as  well.  Figures  22  to  29 
present  the  peak  heating  data  for  boundary  layer  thicknesses  up  to  6 inches. 
At  each  station  downstream  of  the  fin  leading  edge  the  peak  Stanton 
number  correlates  linearly  with  shock  strength  for  a given  boundary  layer 
thickness.  The  slope  of  a line  faired  through  the  data  correlates  with 
"X/6"  as  in  the  peak  pressure  correlations.  This  data  fairing  has  been 

(Stpi/  \ 

Sin  0 j = (0.75,  1 ) in  all  data  sets 

because  it  provides  a simple  yet  accurate  data  fit  for  Sin  0 > 1.3. 

At  small  Sin  0 the  data  actually  approach  (1,1).  The  only  exception 
to  this  data  trend  appears  in  the  Mach  2.95  data  (Reference  14)  in  which 
the  data  tend  to  follow  a line  originating  from  point  (1,1)  even  at 
M^  Sin  0 > 1.3.  No  explanation  for  this  difference  in  data  trends  was 
established.  The  expression  for  the  peak  Stanton  number  derived  by  this 
method  is 


St 

St 


PK 

FP 


(M^  Sin  0 -1)  n^^  + 0.75 


(7) 


where  n^^  is  a function  of  X/6. 


The  function  n^^  is  shown  in  Figure  30. 
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Figure  23.  Peak  Heating  at  Mach  3.00 
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Figure  28.  Peak  Heating  at  Mach  4.75 
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Figure  29.  Pe?.k  Heating  at  Mach  5.04 
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Two  dimensional  interaction  data  has  shown  that  the  peak  Stanton 
number  can  be  correlated  with  the  oblique  shock  pressure  ratio  in  the  form 


“PK  /7(M,Sin  6)2 

Stpp'VP^/  "\  6 } 


(8) 


known  as  the  pressure  interaction  theory  (Reference  15). 


By  plotting  this  Stanton  number  as  a function  Sin  6"  and  calculating 
the  slope  of  the  curve  a value  for  n^^  for  2-D  data  is  derived.  This  value 
is  also  shown  on  Figure  30.  It  can  be  seen  that  within  three  boundary 
layer  thicknesses  of  the  fin  leading  edge  the  3-D  peak  heating  rates  have 
exceeded  the  2-D  values  at  the  same  shock  strength. 

In  the  3-D  interaction  the  region  downstream  of  the  shock  wave  contains 
a vortex  which  dominates  the  flow  field  and  is  not  found  in  the  2-D  inter- 
actions. This  vortex  is  the  phenomenon  which  raises  the  peak  heating 
rates  above  the  2-D  predictions.  The  vortical  amplification  factor  may 
be  found  by  dividing  the  3-D  peak  heating  rate  by  the  oblique  shock  pressure 
ratio  raised  to  the  0.8  power.  This  amplification  factor  was  calculated 
using  the  3-D  peak  heating  correlation  of  Figure  30  and  oblique  shock 
relations. 


Stpk/Stpp  ^ 

(P2/Pi)°’®  ^ 7(M 


n$t  (M^  Sin  9 -1)  + 0.75 


Sin  6)^  -lv°‘® 


(9) 


Equation  9 is  plotted  in  Figure  31  and  shows  the  trend  in  the  amplification 
factor  with  X/6. 


The  location  of  the  peak  pressure  has  been  compared  to  the  location 
of  peak  heating  in  Figure  32  to  see  if  they  are  coincident.  For  X/6  > 10 
the  peak  locations  can  be  taken  as  coincident,  however  in  the  region  at 
small  X/6  the  peak  heating  occurs  slightly  outboard  of  the  peak  pressure 
location. 
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2.  HEAT  TRANSFER  DISTRIBUTIONS 

Token  (Reference  2)  derived  equations  based  upon  stagnation  point 
flow  which  govern  the  heat  transfer  distribution  between  the  peak  and  the 
shock  locations.  The  result  of  his  analysis  was 


, 0.8 
(1  + C J-) 

’*'PK 


(10) 


where  is  the  independent  coordinate  measured  from  the  shock  wave  toward 
the  fin  in  a direction  normal  to  the  fin.  "C"  is  a pressure  gradient 
parameter  defined  by 

0 = "“PK  - (”) 


and  "n"  is  a velocity  gradient  parameter.  In  his  derivation  Token  assumed 
a variation  of  the  velocity  gradient  with  i|>  of  the  form 

6 = K i})*'  (K  constant)  (12) 

The  exponent  in  Equation  12  is  the  velocity  gradient  parameter  in 
Equation  10- 

Through  the  correlation  of  his  data  Token  observed  values  for  these 
parameters  of 

n = 0.5;  0.2  < C < 1 

where  the  weighted  average  for  C was  0.41.  Using  these  values  Token 
successfully  correlated  his  data  taken  at  Mach  3.71  and  Re  ® 3.5  x 10® 
per  foot.  These  values  however  are  not  sufficient  for  the  general 
correlation  of  the  data  presented  in  this  report.  For  large  values  of 
the  shock  strength  the  pressure  gradient  parameter  is  a strong  function 
of  "X/6"  and  "M^  Sin  8"  as  shown  in  Figure  33.  In  this  figure  the 
pressure  gradient  parameter  was  calculated  from  the  data  base  and  plotted 
as  a function  of  M^  Sin  9 . At  a fixed  shock  strength  the  magnitude  of 
"C"  was  observed  to  increase  with  X/6  but  the  precise  function  with 
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X/6  could  not  be  defined.  The  upper  and  lower  curves  on  Figure  33 
approximate  the  maximum  and  minimum  functions  with  X/6  covered  by  the 
data. 


Figure  34  demonstrates  the  effect  of  C on  Equation  10.  It  shows 
that  the  large  values  of  C predicted  by  Figure  33  produce  unrealistic 
heating  rates  for  n = 0.5.  The  value  of  the  velocity  gradient  must  then 
have  a variation  range  similar  to  that  of  the  pressure  gradient.  The 
value  of  n for  any  particular  interaction  must  be  found  by  trial  and 
error  as  shown  in  Figure  35.  The  heating  distribution  is  plotted  according 
to  the  parameters  of  Equation  10,  the  value  of  C is  estimated  from 
Figure  33,  and  values  of  n were  tried  in  Equation  10  until  a reasonable 
fit  was  obtained.  Figure  35  shows  that  n may  be  as  large  as  1.75. 

Equation  10  may  be  used  to  correlate  the  heating  distributions, 
however  to  be  used  as  a prediction  tool  the  shock  heating  rate  must  be 
known.  Correlation  of  the  heating  rate  at  the  shock  location  was 
successful  only  for  M Sin^  0 < 1.7  as  shown  in  Figure  36.  Above  Sin  6 = 

1.7  the  data  trends  diverge  in  a manner  which  could  not  be  accounted  for 
with  the  existing  data  base. 
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SECTION  IV 
CONCLUSIONS 


(1)  It  has  been  shown  that  the  peak  pressure  in  a 3-D  interaction 
can  be  predicted  by 

P nw  0 

#=(M„Sine)P  (I) 


where  Op  is  a function  of  X/6  and  has  a value  of  2.4  at  large  X/6.  Since 
Op  = 2.14  for  oblique  shock  calculations  the  3-D  interaction  can  produce 
peak  pressure  greater  than  oblique  shock  predictions. 


(2)  Once  separation  has  been  established  a plateau  region  develops 
within  the  pressure  profile  which  is  similar  to  that  in  the  2-D  interaction. 
2-D  and  3-0  plateau  pressure  correlate  well  and  separation  occurs  at  73% 

of  the  plateau  pressure  rise.  The  location  of  separation  can  be  predicted 
by  a set  of  hyperbolic  curves  defined  by  Figures  9 to  12. 

(3)  The  peak  heating  rates  in  the  3-D  interaction  are  predicted  by 


^ (M„  Sin  0 -1)  + 0.75 


where  n^^  is  a function  of  X/6.  At  large  X/6  the  value  of  n^^  is  about 
4.75.  Since  2-D  correlations  provide  an  n^^  of  2.65  it  can  be  seen  that 
the  3-D  interaction  produces  heating  rates  which  are  much  more  severe 
than  those  of  a corresponding  2-D  interaction. 
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